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Abstract. Motivated by the central limit problem for convex bodies, we study normal 
approximation of linear functionals of high-dimensional random vectors with various types 
of symmetries. In particular, we obtain results for distributions which are coordinatewise 
symmetric, uniform in a regular simplex, or spherically symmetric. Our proofs are based on 
Stein's method of exchangeable pairs; as far as we know, this approach has not previously 
been used in convex geometry and we give a brief introduction to the classical method. The 
spherically symmetric case is treated by a variation of Stein's method which is adapted for 
continuous symmetries. 



1. Introduction 

Given a random vector X = (X%, . . . , X n ) G M. n , n > 2, and a fixed 9 G S n ~ l , consider the 
random variable 

W e =(X,0), 

where (•,•) is the standard inner product on R". A typical example of interest is when 
X is distributed uniformly in a convex body. In this paper we are interested in determin- 
ing sufficient conditions under which Wg is approximately normal, and in obtaining specific 
error estimates, possibly depending on 9. To do this, we apply Stein's method of exchange- 
able pairs. This technique has not previously been used in studying problems from convex 
geometry, and we believe it will continue to be useful in that context. 

To begin with, we will assume that X is isotropic, that is, that Wg has mean and variance 
1 for every 9 G S' n_1 . Equivalently, X is isotropic if 

EXi = and EX { Xj = 

where 5ij is the Kronecker delta. This is no real restriction, since every random vector with 
finite second moment which is not supported on a proper affine subspace has an afline image 
which is isotropic. 

In the case that the components of X are independent, bounds on the distance of Wg 
from normal follow from classical results. For example, the Berry-Esseen theorem for sums 
of independent, non-identically distributed random variables implies that if X is isotropic 
with independent components, then 

n 

sup \F[Wg <t]- $(t)| < 0.8 f max E^ 3 ) V \9A 3 , 

teR V l<i<n J ' 

2=1 

where $ denotes the standard normal distribution function. There is also a body of work 
going back to Sudakov [3T] (see jH] for a recent contribution and further references) on 
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randomized central limit theorems, which show that under quite general assumptions on X, 
Wg is close to some average distribution for most 9. In these results the average distribution 
may or may not be normal, and "most" may refer to the rotation invariant probability 
measure on S 11 ' 1 or to some other distribution on weights. 

Our motivation in studying this problem comes in part from the so-called central limit 
problem for convex bodies, which is to show that the uniform measure on any high-dimension- 
al convex body has some one- dimensional projection which is approximately normal, or that 
most one-dimensional projections are approximately normal. Most of the results on this 
problem [TJ 12HI prove some form of the latter conjecture (under appropriate assumptions 
on the convex body), and thus fit into the framework of randomized central limit theorems; 
none of the results in these papers identify any specific 9 for which Wg is approximately 
normal. The paper ^U] studies approximate normality of Wg for specific 9 when X is 
uniformly distributed in a cube, Euclidean ball, crosspolytope, or simplex, but in the last 
two cases only for a very restricted set of 9 and with rather limited quantitative information. 

Of course, there is no hope to identify any specific 9 for which Wg is approximately normal 
without some additional assumptions on the distribution of X. Here the additional assump- 
tions we consider are more geometric than probabilistic in nature. Specifically, we consider 
distributions which have a sufficiently rich class of symmetries, although we emphasize that 
our results do not require X to be drawn from a convex body, or even to be continuous. 
Stein's method, described in Section El allows us to take advantage of these symmetries in 
order to reduce normal approximation to estimation of certain low-order moments. 

Our first main result treats distributions which are symmetric with respect to reflection in a 
suitable collection of hyperplanes. Our hypothesis encompasses both the class of distributions 
which are coordinatewise symmetric ( Corollary |2J) and those with the symmetries of a regular 
simplex (Corollary The error bounds for the approximations are in many cases small 
enough to derive multivariate randomized versions which improve on existing results. 

Our second main result, treating spherically symmetric distributions, is proved by a vari- 
ation of the classical version of Stein's method adapted to take advantage of continuous 
symmetries. This result has as corollaries several classical results as well as some new ap- 
plications. We also make connections with Poincare inequalities for probability measures on 



The layout of this paper is as follows. We first define notations which will be used through- 
out the paper. In Section |21 we state our results and several corollaries, and give comparisons 
to existing results. Section El gives a brief introduction to Stein's method. Section 0] contains 
the proofs of the first main theorem and its corollaries. Section El contains the proof of the 
second main theorem using the variation of Stein's method described above, and the proofs 
of its corollaries. 



Notation. Let £™ = (W 1 , || • \\ p ), where || • || p denotes the norm 
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for 1 < p < oo, and 

l<i<n 

For v G R n , define v ® v : R n -> M n by 

t> £g> v(x) = (x, v) v; 

if v G S" 1-1 , u ®t> is the orthogonal projection onto the span of v. A set of vectors u\, . . . , u m G 
S 11 - 1 such that 

m 

(i) y Ui®Ui = —i n , 

is known in the signal processing literature as a normalized tight frame (J„ is the identity 
on R n ). By taking the trace of both sides of (JTJ), one can see that — is the only possible 
constant that can appear. 

The Grassmann manifold of /c-dimensional subspaces of M. n is denoted G n ^\ it is equipped 
with a unique rotation-invariant probability measure X n ^- For a fixed subspace E C M n , let 
Pe denote the orthogonal projection onto E, and 7^ the standard Gaussian measure on E. 

For a random variable or random vector X, let £>(X) denote the distribution of X Given 
two probability measures \i and v on E, define the T-distance between them as 

T(fj,,u) = sup{J/x(iZ) — v{H)\ : H is an affine half-space of E}. 

In particular, 

T(L(P E (X)), lE ) = sup \¥[W e < t] - $(t)\. 

eeEns 71 - 1 , teR 

This is a geometrically natural notion of distance between measures since it is invariant 
under nonsingular affine transformations and is thus not tied to any particular coordinate 
system. In addition, the topology induced by T on the space of probability measures on E 
is stronger than the w* topology. 

The total variation distance between two probability measures /jl and v is 

dxvit 1 , 1 ') — 2 sup{ |//(^4) — u(A)\ : A is measurable}. 

Recall that if (jl and v both have densities, then their total variation distance is the Li 
distance between their densities. 

Finally, note that symbols like c, c\, etc. which represent absolute constants may have 
different values from one appearance to the next. 

2. Statements of results 

Theorem Q the first main result of this paper, is based on existing normal approximation 
results proved via Stein's method. Corollaries HJ 01 EJ and HI are all applications of Theorem 

12 

Theorem 1. Let Ui, . . . ,u m G >S n ~ 1 be a normalized tight frame, and for any x G M. n let 

X(i) = (x, Ui) . Suppose that X is a random vector whose distribution is invariant under 
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reflections in each of the hyperplanes uf. Let 6 G S 71 ^ 1 be fixed. Then 



sup|P[W e <t}- $(t)| < 2 

teR 



\ 



n 2 
m 2 



//4 

7T 



//in addition max |X(j)| < a almost surely, then 

Ki<m 



suplPfWfl < t] - < 24 



\ 



^ .EpC 2 , X?. 1 -1 



+ 172na 3 max |#m| 3 . 



Ki<m 



The constants that appear in these statements, and explicit constants which appear in 
any of the results that follow, are generally not the best possible, and are included only for 
concreteness. 

It is not obvious from the statement that the error estimate in Theorem is useful. 
However, the proofs of Corollaries 121 and El will show that Theorem ^ allows different cases 
of geometric interest to be treated easily in this unified framework. 

Borrowing terminology from the geometry of Banach spaces, call X unconditional if its 
distribution is invariant under reflections in the coordinate hyperplanes, or equivalently, if 
X has the same distribution as {s\Xi, . . . , e n X n ) for any choice of e G {—1, 1}™. 

Corollary 2. Let X be unconditional and isotropic, and let 9 G S"™ -1 be fixed. Then 



suplPfWe < t] -$(t)| < 2j( maxEX 4 4 )||fl||4 + maxCov(X 2 ,X 2 ) 

V V l<i<n J i^j ' ■> 



1/4 , 

I max 

7T / \l<*<n 



If moreover X G [—a, a]™ almost surely, then 



sup|P[Wfl < t] -$(t)\ < 2AM max EXf) \\6\\\ + max Cov(X?,Xf 
teR " V ^ ' &i 

+ 172na 3 ||fl| |3 



I oo ' 



The statement of Corollary ^motivates the introduction of the following definition, taken 
from [23] ■ A random vector X has the square negative correlation property if 

EX 2 * 2 < (EX 2 ) (EX 2 ) forz^j, 

i.e. if Cov(X 2 , X 2 ) < for % ^ j. 
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A random vector X is called log-concave if it has a logarithmically concave density with 
respect to Lebesgue measure on IR n . In particular, a random vector which is uniformly 
distributed on a convex body is log-concave. 

It natural to conjecture (cf. Section 5 of [23]) that an isotropic unconditional log-concave 
random vector has the square negative correlation property. However, this is not the case, 
as shown by the example [3] of the density 

(2) a n e~ bnM ^ 

on R n , where a n and b n are appropriate normalizing constants. (Counterexamples also exist 
which are uniformly distributed in a convex body 0.) However, the weaker conjecture that 
under these conditions 

(3) Cov(X?,X])<- for i ^ j, 

J n 

for some absolute constant c, is open (and is in particular satisfied by X with the density 
(0)). This conjecture is related to the Kannan-Lovasz-Simonovits conjecture on isoperimetric 
constants [T7] (see [S] for a discussion of this issue, and cf. Corollary [TT] below) . 

The error bounds in Theorem ^ are small enough in many cases to show that Wg is 
uniformly close to normal for all unit vectors 9 in a typical subspace E C M. n of relatively 
large dimension. This means that the projection of X onto E is close to normally distributed 
in the sense of T-distance. In order to quantify this phenomenon, for e > 0, define 

A ntk (e) = {Ee G n>k : T(tL(P E (X)), lE ) < e}; 

that is, A n>k (e) is the set of ^-dimensional subspaces E C lR n such that the projection of X 
onto E is e-close to normal in the sense of T-distance. 

The following lemma, proved in allows normal approximation in the sense of total 
variation to be deduced from approximation of distribution functions in the case of log- 
concave distributions. 

Lemma 3 (Brehm-Hinow-Vogt-Voigt). There is an increasing function j3 : (0,oo) — ► (0,2] 
satisfying 

(5{t) = O (y~t\ogtj ast^O 
such that for any log-concave random variable W, 

drv(HW),-yn) < [3 ( sup \F[W < t] - $(t)\) . 

\tm / 

All references to (3 in the next three results are to the function in Lemma El 

Corollary 4. Let X be unconditional and isotropic, with the square negative correlation 
property. Then there are constants ci, . . . ,Cn, independent of X and n, such that each of 
the following holds. 

(1) IfE\Xi\ 3 < a and EXf < b for all i, then 

sup\F[W e <t)- $(i)| < 2(Vb\\9\\l + v^lla 72 ) 

teK 
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for all 9 G S 12 ^ 1 , and 

X n , k (A n>k (e)) > 1 - Cl exp [-c 2 min{cr 2/3 £ 4/3 , b- 1/2 e}n] 
for e > c 3 max{^/an~ 1 / 4: 1 \/bn~ l l 2 } and k < c 4 min{a _2//3 £ 4 / 3 , b~ x l 2 e\n. 



(2) //X zs log-concave, then 



sup|P[jy e < t] - < c 5 \\9 



1 3/2 



/or a// G S"^ 1 , and 

An,fc(Ai,fc(e)) > 1 - ciexp [-c 6 e 4/3 n] 
for e > c-jn~ x l^ and k < Cg£ 4//3 n. Furthermore, 

d T v{£(W e ), lR ) < (3(c 5 \\ef 3 /2 ) 

for all Be S n -\ 
(3) If X G [— a, a] n almost surely, then 

sup|P[W fl <t]- < i96na 3 ||0||^ 

*eR 

/or a// G S 1 " -1 , and 

X n ,k{A n ,k{ £ )) > 1 - c i ex P [-c 9 a~ 2 e 2/3 n 1/3 ] 

fore > Cioa 3 (logn) 3 / 2 n _1//2 and k < cua~ 2 e 2 ^n 1 ^. If moreover X is log-concave, 
then 

«W(£(V^), 7 r) </?(l96na 3 ||0||L) 

/or all 6 G S"- 1 . 



The square negative correlation property is included as a hypothesis of Corollary 0] only 
for convenience. Replacing it with the hypothesis Q would result only in a weakening of 
the constants that appear, and in fact the even weaker hypothesis 

Cov(X 2 ,X 2 )<^ fori^j 

would suffice for the same conclusion in part|2l In particular, the conclusion of part |2] applies 
to X distributed according to the density (J2J). 

One could also deduce randomized total variation results for one- dimensional projections 
in parts 121 and El of Corollary 0J but the statements are more complicated. 

Naor and Romik [23 Theorem 5] proved a result comparable to the randomized statement 
in part ^ of Corollary HJ Under similar hypotheses (but without unconditionality), they 
showed 

\ n ,k{A n ,k{£)) > 1 - ^-exp [-c 2 b' 1 e 4 n] 

for e > and k < csb~~ 1 e 4 n, so Corollary 0] improves on the dependence on both b and e in 
the unconditional case. In the case that X is uniform in a convex body (hence log-concave) 



THE CENTRAL LIMIT PROBLEM FOR RANDOM VECTORS WITH SYMMETRIES 



7 



and has the square negative correlation property, Antilla, Ball, and Perissinaki jH Theorem 
2] showed 

(4) X n ,i(Ai,i( £ )) > 1 - nexp [-ae 2 n] 

for e > C2 r n~ 1 / 3 . Part 121 of Corollary 0] improves (in the unconditional case) on their depen- 
dence on e, although for a slightly more restricted range of e, and does not require that X 
be chosen from a convex body. In the case of certain bounded distributions, part H3 improves 
further on the results of [11123, as wm be illustrated by Corollary El below. 

The next corollary treats a class of examples of particular interest in asymptotic convex 
geometry, namely, X chosen from various natural distributions on the unit balls of the spaces 
£p. In addition to the uniform measure on the interior, there are two geometrically natural 
measures on the boundary of a convex body K whose interior contains the origin. First 
there is (n — l)-dimensional Hausdorff measure, or surface measure. Second, there is cone 
measure /x, defined by 

H(A) = vol ( |J tA) for A C dK. 

*6[0,1] 

Cone measure is the measure on dK for which there is a straightforward extension of the 
familiar polar integration formula, with dK replacing S 11 ^ 1 . 

Corollary 5. For 1 < p < oo, let X have one of the following distributions: 

(1) uniform measure on the ball of £™, scaled to be isotropic; 

(2) normalized cone measure on the sphere of £™, scaled to be isotropic; or 

(3) normalized surface measure on the sphere of £ p , scaled such that the normalized cone 
measure is isotropic. 

Then there are absolute constants c\, . . . , c 5 and constants di )P , . . . , d^ p depending only on p 
such that 

(5) sup|P[W e <t]- < min UpfJ 2 , d 1}P n 1+ l \\9\\ 3 \ 

teR *> J 

for all Be S n ~ l , 

X n> k(Ai,k( £ )) > 1 ~ c 2 exp [-c 3 e 4/3 n] 
for e > CiU~ l / A and k < c^e^n, and 

K,k{ A n,k{ e )) > 1 - c 2 exp 

for e > d 3 ^(\ogn) 3 / 2 n~p~2 and k < d^e 2 ! 3 ^'? . 

Furthermore, in the case that X is chosen uniformly from the reseated £™ ball, 

drv(^(We),^) < (3 (min {dpfj 2 , d hp n 1+ l ||0||^}) 

for allOe S n -\ 

To compare the two bounds in (j5J, note that since ||0||oo < ||#||3 < W^IU = 1 and ||#||3 > 
n -1 / 6 , the first error bound is better for all 9 when 1 < p < 4 (ignoring the constant factors). 
On the other hand, for the principal diagonal 9 = n~ x l 2 ^ e« (which roughly captures typical 



-d 2 , p e 1 n 3 p 



8 



E. MECKES AND M. MECKES 



behavior for norms on S n x ) the second bound is better for p > 12. In particular, for 
p > 18, Corollary El improves on the typical rate of convergence to normality of about n" 1 ' 3 
which follows from . Corollary also improves on Theorems 7 and 8 of , which show 

A„ )fc (A„ )fc (£)) > 1 - ^ exp [-c 2 e 4 n] 

for e > and k < c^e A n in cases El and El of Corollary 

Brehm and Voigt JU] considered X uniformly distributed in the rescaled £™ ball for p = 
1, 2, oo. See the discussion of Corollary |U] below for the case p = 2. In the case p = oo, they 
derive sharper error bounds for general 9 than here. In the case p = 1 (the crosspolytope) 
however, they consider only the case 9 = n' 1 ^ 2 ^ Ci for n — > oo, and do not obtain an explicit 
rate of convergence, so Corollary provides a substantial generalization and strengthening. 

As discussed earlier, the form of Theorem ^ is general enough to accommodate the sym- 
metries both of unconditional distributions and of a regular simplex, which is treated in the 
next result. 



Corollary 6. Let A n = ^Jn(n + 2) conv{i>i, . . . , v n+ {\ be an isotropic regular simplex, where 
Vi G S 11 ^ 1 . Let X be uniformly distributed in A n , and let 9 e S 1 ™ -1 be fixed. Then there are 
constants ci, . . . , c 5; independent of n, such that 



sup|P[W fl <t\- < ci, 



n+1 



\ i=i 



and 



n+1 



\i i=i 

/or ai/de and 

\ nt k(A ntk {e)) > 1 - c 2 exp [-c 3 e 4/3 n] 
/or e > c^n" 1 / 4 and k < c^e^n. 



The case in which X is uniformly distributed in a regular simplex was also considered 
in ^Uj. However, the results there consider only a certain discrete set of 9 (roughly those 
for which the behavior of Wg is best) for n — > oo, and do not derive any explicit rate of 
convergence to normality. 

The remaining results are not based on existing normal approximation theorems; instead, 
the proofs use a variation of Stein's method of exchangeable pairs, adapted to situations in 
which there are continuous symmetries. This variation was introduced by Stein in [30J and 
developed further by the first-named author in [21] in studying functions on the classical 
matrix groups. 
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Theorem 7. Let X be an isotropic random vector, with finite third moment, whose distri- 
bution is spherically symmetric. Then for any 9 G S 1 " -1 , 

<M£(W fl ),7te) < 4E|1 -E[X 2 2 |X! 



4 9 

< E X 2 — n\ + — 



<-^ryJv°*(\\x\\*) + - 

n — 1 v n 



The latter bounds in Theorem [7| reduce normal approximation of one-dimensional projec- 
tions of spherically symmetric random vectors to the problem of estimating deviations of 
|| X HI from its mean. Some kind of such concentration of ||X|| 2 is either explicitly a hypoth- 
esis, or closely related to the key hypothesis, in many of the existing results on the central 
limit problem for convex bodies, cf . [U IE1 HE HE] • 

Before proceeding to some classical consequences of Theorem we first state a corol- 
lary giving a new randomized central limit theorem, which in particular gives information 
about the central limit problem for convex bodies. Suppose X is a (not necessarily spher- 
ically symmetric) isotropic random vector and U is a random n x n orthogonal matrix, 
distributed according to Haar measure independently of X. Notice that the distribution of 
W = (UX, ei) = (X, U" 1 ei) is the average (with respect to the rotation invariant probability 
on S 71 ^ 1 ) of the distributions of Wg over all 9 G S n ~ l . As mentioned in the introduction, the 
distribution of W is the object of some of the work on randomized central limit theorems. 
Since X = UX is a spherically symmetric isotropic random vector and ||X||2 = ||X||2, the 
following is an immediate consequence of Theorem [7| 

Corollary 8. Let X be an isotropic random vector with finite third moment, and let W be 
as defined above. Then 

d Tv mW),<n) < -^-E\\\Xg - n\ + 

n — 1 n 



<^-r\fv™(\\X\\t) + ^-r 
n — 1 v Ti — 1 



In the case in which X is uniformly distributed in if C 1", the density of W gives the 
average (n — l)-dimensional volume of a hyperplane section of K at a given distance from 
the origin. For K a convex body, Bobkov and Koldobsky jH] proved a pointwise bound on 
the difference between the density of W and the standard normal density, with a bound 
which also explicitly involves the variance which appears in Corollary |S] (and which is also 
essentially of the order n -1 as long as the variance is not too big). Corollary |H] gives instead 
an L\ bound on the difference of these densities. See ^H] for an earlier asymptotic result 
in the case that K is a cube, and jl] for a generalization of the result of :8j to arbitrary 
distributions and a multivariate version for sections by /c-codimensional affine subspaces. 

The following easy corollary of Theorem [7| is well-known; versions for higher-dimensional 
projections are proved in ^S] for the sphere and in ^D] for the ball. Theorem [7| allows the 
cases of both the ball and the sphere to be presented simply as part of a unified framework. 
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Corollary 9. If X has the uniform distribution on the Euclidean ball of radius \Jn + 2 or 
the uniform distribution on the sphere of radius \/n, then for any 9 G <S n_1 

a 

d TV (H(W e ),j R ) < -, 

n — 1 

where a = 16 in the case of the ball and a = 8 in the case of the ball. 

Corollary El gives the correct order of approximation in both cases, although the constants 
can be improved. 

The first error estimate in Theorem [7| is also strong enough to recover, as an immedi- 
ate consequence, a version of the characterization of the normal distribution as the unique 
spherically symmetric product measure on R n . (The two-dimensional version of this charac- 
terization is the classical Herschel-Maxwell theorem; see ^T] for various other versions.) 

Corollary 10. A random vector X with finite third moment has the standard normal dis- 
tribution if and only if X is spherically symmetric and has independent components with 
variance 1. 

Recall that a random vector X is said to satisfy a Poincare inequality with constant Ai 
(the spectral gap of X) if 

(6) A x Var(/(X)) <E||V/(X)||^ 

for every smooth / : W 1 — > M.. The last estimate in Theorem [7| provides a connection between 
normal approximation and spectral gap estimates. A similar connection has been observed 
in a different but related context by Bobkov and Koldobsky [S]. 

Corollary 11. Let X be an isotropic spherically symmetric random vector with spectral gap 
\\. Then for any 9 G i?" -1 , 

drv(£(Wfl),7te) < -j==. 

One concrete application of Corollary ITT1 is the following. 

Corollary 12. Let X have the isotropic spherically symmetric exponential density 

n P -MM| 2 

where a n and b n are appropriate normalization constants. Then 



d T y(£(iy e ),7 M ) < 



n 



1/2 



Bobkov [7j showed that for this distribution, ^ < Ai < 1, and so Corollary^Jis immediate 
from Corollary^! Since this distribution is given explicitly, one can obtain an error estimate 
of the same order by directly estimating the variance in Theorem however, there is a large 
literature on spectral gap estimates in much less explicit contexts using only certain geometric 
assumptions, typically diameter and/or curvature bounds (see [12] for a survey and further 
references). Corollary ^2 thus allows the treatment of distributions about which one has 
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geometric information resulting in spectral gap estimates, but for which direct computation 
of the variance term in Theorem [7| is not possible. 

There is also the following complex analogue of Theorem All of the previously defined 
notation of this paper used here should be reinterpreted for vectors in C n in the most obvious 
way. 

Theorem 13. Let X G C n be a random vector with finite third moment such that EX; = 



for each i, EXiXj = EXiXj = if i ^ j, and E(Re X, 



E(ImX; 



1 . Suppose 



the distribution of X is invariant under multiplication by a unitary matrix. Then for any 
9eS, 



n-l 



d TV {L{ReW e ),^) < 4E 



< 



1 

2n 



n 



2(n+l) 



E[|X 2 | 2 |Xj 



n 2 — 1 



Var 



1*111) + 



n 



n 



Note that under the hypotheses of Theorem E3J E[|X 2 | 2 |X!] = 2E[(Re X 2 ) 2 \X 1 ], and so 
the appearance of the factor of 2 inside the first bound above is to be expected. 



3. Background on Stein's method 

The essential idea of Stein's method is the notion of a characterizing operator. Say that 
T is a characterizing operator for a distribution on R if the following conditions hold: 

(1) J T f{t) dfi(t) = for all / such that T f is /x-integrable, and 

(2) if v is a probability measure on R such that J T f{t) dv{t) = for all / with T Q f 
z^-integrable, then /j, = v. 

A characterizing operator is a strong characterization of a distribution, in the sense that if 
T Q is characterizing for \i and v is a measure such that J T Q f(t) du(t) is small for a large 
class of test functions /, then u « fi in some sense. 

This idea is quantified in the method of exchangeable pairs as follows. Let W = W(u) be 
a random variable defined on a probability space (Q, P), and let X be a space of measurable 
functions on Q. Think of E as a linear map from X to IR, E/ = J f(u>) dP(u). Let X Q be 
a space of measurable functions on R and let E D be the linear function on X Q defined by 
E / = / f(t) dfj,(t) for some fixed measure \i. In our applications, fi will be the standard 
normal distribution, but one of the advantages of Stein's method is that the set-up is quite 
general and can be adapted to various other measures. The random variable W induces a 
map (3 : X Q — > X defined by 

(7) (3f(u) = f(W(u)). 

Now, construct a symmetric probability Q on Q x Q with margins P (i.e., Q(A x B) = 
Q(B x A) and Q(A x fi) = P(A)). Note that this is the same as constructing an exchangeable 
pair (W, W) = (W(u), W(u/')) from W. Let J be a space of measurable, antisymmetric 
functions on Q x Q and use Q to define a map T : "J — > X by 

(8) Tf{u) = E Q [f(u;,u')\Lo]; 
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note that by exchangeability and antisymmetry, ET = on 3 r . Let be another space of 
measurable functions on R, possibly the same as X Q , and let T a : J Q — > X Q be a characterizing 
operator of /i. Let U Q : X G — > Jo be a pseudo- inverse to T Q , in the sense that 

(9) T UJ(t) = fit) - EJ. 

Finally, let a : £F — > 3 r . All of these definitions are summarized in the following diagram: 
(10) 




The following easy lemma (SHj is the quantitative version of the heuristic at the beginning 
of the section. 

Lemma 14 (Stein). Suppose that in the diagram of spaces and maps above, ET = and 
T U = Id-E G . Then 

E(3-E =~E((3T -Ta)U . 



Note that for / e X Q , 

Epf-EJ = Ef(W)-Ef{Z), 
where Z is a random variable with distribution \l. The strategy is apply Stein's lemma to 
bound this difference uniformly over a large class of test functions. 

To apply the method of exchangeable pairs, one needs an approximating distribution and 
a characterizing operator. In this paper, the measure \x will be standard Gaussian, with 
characterizing operator 

TJ(t) = f(t)-tf(t) 

and pseudo-inverse 



u f(t) 



1+2 

e2 



* 2 ! [f{s)-E f]e-^ ds. 

J — oo 



That E T = can be verified by integration by parts, and verifying T U = Id— E Q is just 
calculus. 

The following estimates for U Q are proved in P- 25] and are useful in estimating the 
error term from Lemma ITU 



(11) 

(12) 
(13) 



\U OJ 



<J-\\f-E f\\ 00 <V2n 



||(f/o/) , ||oo<2||/-E /|| 00 <4 

iKt/o/nu^u/'iu 



Using this general set-up one can prove the following abstract normal approximation 
theorem, which is the main tool used in the proof of Theorem The first statement was 
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proved by Stein in j2Hl Lecture III]; the second statement was proved by Rinott and Rotar 
in EE]. 

Proposition 15 (Stein; Rinott-Rotar). Let (W, W) be an exchangeable pair of random 
variables such that 

EW = 0, EW 2 = 1, 

and 

E[W-W'\W] = XW 

for some A £ (0, 1). Then 

sup|P[W<t]-$(t)| < T \/Wp - + (2ny 1/ \ -E\W - W'\ 3 . 

tm A V A 

If moreover \ W — W'\ is almost surely bounded, then 

12 43 
sup|P[iy < t) - $(i)| < — ^/Y^E^T^Wf\W] + — \\W- w'wic. 

tern. A A 

Remarks 

(1) For Stein's lemma to hold, no assumptions are needed on the map a : 3 Q — > 3 other 
than that it does in fact produce antisymmetric functions. For normal approximation, 
one usually uses the map 

af(u,u/) = a(W{J) - W{u)) [f(W(co')) + f(W(u))}, 

for some suitable choice of a. 

(2) In the description of the method, we have been quite cavalier about exactly which 
function spaces should be used. Of course, the choice of the space of test functions 
X determines the type of convergence; in practice, one generally fixes X Q first and 
chooses the remaining spaces in some way which guarantees that all of the maps do 
fit into the diagram as shown. The diagram and Stein's lemma are described here 
mainly as motivation for the approach taken in Section El and are often not used 
directly but guide for how to proceed. 

(3) Proposition has been applied in many different contexts. In particular, Holmes 
and Reinert [13] use it in analyzing statistics similar to our Wg while studying process 
approximation for the bootstrap. 

(4) This paper is the first that we know of to apply the exchangeable pairs approach to 
Stein's method in convex geometry. However, Reitzner [23] used a rather different 
approach to Stein's method (based on dependency graphs) in proving central limit 
theorems for random polytopes. 

(5) We have described here how to apply Stein's lemma to normal approximation; appli- 
cation to Poisson approximation has also been extensively developed, see |12j . 



14 E. MECKES AND M. MECKES 

4. Proof of Theorem [T] and its consequences 

Proof of TheoremUi Let X be a random vector invariant under reflections in the hyperplanes 
defined by a normalized tight frame Ui, . . . , u m . Define an exchangeable pair (W, W) as in 
Proposition El as follows. Let I be chosen uniformly from {1, . . . , m}, independently of X, 
and define 

X' = X - 2X {I)Ul , 

i.e., X' is obtained from X by reflection in the hyperplane uf. Then (X, X') is an exchange- 
able pair of random vectors by assumption. Define W = Wg = (X,8) and W = (X',9). 
Now EW = and EW 2 = 1 since X is isotropic, and 



E[W — W'\W] = E 



W 



1e 

m 



m 

^2ui®Ui(X),9 



4=1 



w 



-w 



n 



by equation (JTJ). 

To apply Proposition ITB1 (with A = -), it remains to estimate the quantities 



V&tE[(W - W') 2 \W], E\W-W'\ 3 , and \\W - W 
(the last in the case that max \ X^\ < a almost surely). First, 

E(E[(W - W') 2 |Wj) = E(E[VF 2 + {W'f - 2WW|W]) = 
and by the conditional form of Jensen's inequality, 



4 

> 

n 



E(E[(W - W') 2 \W}) 2 < E{E[{W - W') 2 \X}) 

= E(E[(2X (I) e (I) ) 2 \X]) 
16 



Yl e h e (i) E [ x h x u)\> 



so 



(14) 



VarE[(W - W') 2 \W] = E(E[(W - W') 2 \W}) 

16 / n 2 



16 

~2 



n 



< 
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Next, 



E\W- Wf = 8E|X (J) (J) | 3 



(15) rn 



X>(,)| 3 E|X ( . 



1 3 



i=l 

m 

|3 



8 / \ 
< - max E|X (J) | 3 V|%)| 

777. V Ki<m w / w 



Finally, if max \Xu\\ < a almost surely, then 

(16) \W -W'\ = 2|0 (/) ||X (/) | < 2a max |%) | 

almost surely. Inserting the estimates (|14|) . (|15j) . and (|16|) into Proposition [TBI now proves 
Theorem ^ D 

Remark: One can also prove, by the same method, a version of Theorem Q for random 
vectors which are invariant under reflections in subspaces of arbitrary dimensions. Rather 
than a normalized tight frame, one would consider a decomposition 

m 
i=l 

for subspaces Ex, ... , E m C M n . In applications this can lead to a nontrivial approximation 
result if the maximum dimension of the remains bounded (or grows very slowly) as n 
grows. This context covers, for example, X uniformly distributed in the £ p -sum of m copies 
of a fixed symmetric convex body K <zM. k . 

Proof of Corollary^ The corollary follows easily from Theorem Q by taking as the normal- 
ized tight frame the standard basis of lR ra (so that xu) = Xi). The form of the second error 
term in both statements follows immediately from Theorem ^ To estimate the first error 
term, note that 

YJ < max E.Y* ^ + E ^ 

i,j = l V ' i=l V ^ ' i,j = l 



m 

max EX. 4 V^ 4 + maxEX 2 X, 2 



l<j<m J I ' * i^j 



since \\6\U = 1, and 



EX?X] = Cov(X, 2 ,Xj) +1 
since X is isotropic. □ 

The following lemma is used to prove the randomized statements in Corollaries 01 
and El It follows easily from a concentration inequality implicitly proved by Gordon 
(see Theorem 6 in (221 f° r an explicit statement), together with the well-known asymptotic 
orders of the averages of PI norms over S 1 " -1 . 
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Lemma 16. There are absolute constants ci, . . . , C4 such that the following hold. 

(1) If 5 > C\n~ x l % and k < C25 2 n, then 

X n , k ({E E G nyk : ||0|| 3 < 5 V0 E E n S n - 1 }) > 1 - c 3 e- c ^ n . 

(2) If 8 > ci??," 1 / 4 and k < c 2 S 2 n, then 

Xn,k({E e G n>k : ||0|| 4 < 6 V0 E E n S n - 1 }) > 1 - c 3 e- C4<52n . 



(3) IfS> Cl J l ^andk< c 2 5 2 n, th 



en 



n— 1\\ \ 1 „ — C4<5 2 n 



A„, fc ({£ G G n , fc : H0IU < 5 V# E E n > 1 - c 3 e 

Proof of Corollary^ 

(1) The first statement is immediate from the first statement of Corollary 121 the second 
follows from parts Q and El of Lemma E3 

(2) It is a well-known consequence of Borell's lemma (cf. [213 Section 2.2]) that there is 
an absolute constant c, independent of X, n, and p, such that 

(17) (E\(X,y)\f) 1/p <cp(E\(X,y)\ 2 ) 1/2 

for any log-concave random vector X, p > 2, and fixed vector y. Thus part applies 
with some absolute constants a and b. Furthermore, since 1 1 6* 1 1 4 < 1 1 6> 1 1 3 < ||0|| 2 = 1, 
the first term in the r.h.s. of part Q is not of larger order than the second term. The 
total variation bound follows from Lemma El and the fact that any projection of a 
log-concave measure is again log-concave [2I|. 

(3) From the second statement of Corollary |21 and the trivial estimate KXf < a 4 we 
obtain 

sup|P[W«, < t] - $(t)| < 24a 2 ||0||^ + 172na 3 ||0||^. 



<n 1 / 4 



The estimates 

l|0||4 ^ >> |"||x 

and 

1 = ||0|| 2 < V^I|0|U 

are well-known, and a > 1 since X is isotropic. Therefore 

a 2 ^||4<a 3 v^l^llL<^ll^l| 3 oo, 

which proves the first statement. The randomized statement follows from part El of 
Lemma El and the total variation bound follows from Lemma El 

□ 
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Proof of Corollary^ The square negative correlation property was proved for the uniform 
measure on the ball of £™ in |3. (see also PQ) and for the cone measure in (221- The uniform 
measure on the ball is log-concave by the Brunn- Minkowski theorem, and it is not hard 
to show (cf. j2Hj) that E|Xj| 4 < c for some absolute constant in the case of cone measure. 
Finally, it is well-known that if the £™ ball is scaled so that its uniform measure is isotropic, 
then it is contained in [ 1//p ]™, where Clp IS cL constant depending only on p; it is 

not hard to show the same is true of the normalized cone measure. 

Using all these facts, the statements for uniform measure on the ball and cone measure 
on the sphere follow from Corollary the total variation bound for the ball follows from 
Lemma E3 

The statements for the surface measure then follow from the fact, proved in [2HJI22] that 
the total variation distance between the cone and surface measures is at most -§= for some 

s/n 

absolute constant c, and both of the error estimates are of at least this order (cf. the proofs 
of El and H of Corollary ij) . 

□ 



Proof of Corollary \Q First, if A n = ^Jn(n + 2) conv{t>i, . . . , v n+ i} is a regular simplex, then 
the Vi form a normalized tight frame and also satisfy 

n+l 

(18) 5> = 

i=i 

Both of these facts are well-known and can be seen as consequences of John's theorem on 
contact points between a convex body and the minimal volume ellipsoid containing it [To] . 
To see Corollary El as a consequence of Theorem ^ consider the vectors 



/ Tl 

It is not hard to show from (|18|) that ||wy||2 — 1 for each i ^ j, and that the Uij form a 
normalized tight frame because the Vi do. Reflection in is a reflection which interchanges 
the vertices in the directions Vi and Vj and leaves all other vertices of A n fixed. Theorem 
[^can thus be applied (with m = n(n + 1)), provided that X is indeed isotropic under this 
scaling. This is not obvious at this point, but follows easily from ()19|) below. 

In order to compute the relevant expectations, one can embed A n isometrically in M n+1 
by the affine map with 

y/n{n + 2)vi i — ► y/(n + l)(n + 2)ei; 
the image of A n under this map is 

A' n = y/(n+l)(n + 2) conv{ei, . . . , e n+1 }. 
Let Y be the image of X under this isometry; Y is uniformly distributed in A' . Then 
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and so, adapting the notation of Theorem [TJ 

X(ij) = (X, Uij) 



n 



2(n+l) 



7? 



((X,v t )-(X lVj )) 



2(77+1) 



n 



Y,J n -±±e, 
n 



V2 



{Yi-YiY 



The joint moments of the Yi are given by 

"n+l 



(19) 



E 



i=i 



[(77 + l)(77 + 2)] r/2 _r7!^ T 1 



(n + r)! 



i=i 



where r = ^2 r j . This formula follows easily from Lemma 1 of [21] • Using this one can show 
that EX ( 2 . j} = 1 and 
i 7^ j and k ^ I, then 



that HL = 1 and thus that the stated normalization for A n is correct. In addition, if 



(20) 

and 
(21) 



Fr2 Y 2 _ (n + l)(n + 2) 



(77 + 3)(77 + 4) 



i if {£,i}n{fe, = 0, 

3 if |{i,j}n{fc,Z}| = l, 

6 if \{i,j}n{k,i}\ = 2, 



EIX 



< 3^ 



V(n+l)(77 + 2) 



77 + 3 



< 3^. 



(What is really needed about this latter quantity is just that it is bounded by an absolute 
constant; this also follows immediately from Borell's lemma Q17]).) 
To estimate the first error term from Theorem^ by (|20|1 . 



Yl 6 h) d (ki) E [ X (ij) X (ki)] 
(«),(fcO 

(77+ 1)(t7 + 2) 



(t7 + 3)(t7 + 4) 



4 



(y) 



_(y),(fci) |{i,i}n{fc,«}|=i 

where in all of the above sums, indices run from 1 to n + 1, % ^ j, and k ^ I. Since the in- 
form a normalized tight frame, 

(y) 

and by the Cauchy-Schwarz inequality, 



|{i,j}n{fc,z}|=i 



# 2 fl 2 <-" 



|{»j}n{fe,J}|=i 
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By ()18p and the fact that the Vi form a normalized tight frame, 

/ \ 2 n+1 

E4-)=(-^) E ««>«<>-<«. •*»' 



2 n +! 



2n+_ 
v ' i=i 



and so 
(22) 



^ E W[4)4o]- 1 S''D«.» i > 4 + (^ 

; (ii).(kl) i=l V 7 



(«),(«) 

To estimate the second error term from Theorem ^ 

1/3 , /n+1 



(23) 



3 



< 



by the triangle inequality for the norm. The first statement of the Corollary now 

follows by inserting (|2ip. ()22|). and (|23p into Theorem ^ and noting as in the proof of part 
121 of Corollary 0] that the first error term is of smaller order than the second error term. The 
total variation bound then follows from Lemma E3 

The randomized statement essentially follows from Lemma HKl as in the proof of Corollary 
I3J In this case one actually needs not part Q of Lemma EE but the same estimate for the 
norm 

(n+1 \ X /3 

i=i 

which can be proved in the same way. □ 



Remarks: The same issue of estimating covariances of the squares of frame components 
of X, which is explicit in the error bound of Corollary |21 also arises implicitly in the proof of 
Corollary IH1 In the latter case Cov(X£~, is not always small; the key point is that it is 

at worst a positive constant, and this only happens for a negligible fraction of pairs (ij), (kl). 
Also, it is clear that the proof of Corollary |H1 can be adapted to treat other distributions 
which posses the same symmetries as a centered regular simplex. 

5. Infinitesimal rotations 

This section is mainly devoted to the proof of Theorem [7| In the usual exchangeable 
pairs approach to Stein's method described in Section one starts with a random variable 
W and makes a small change to get a new random variable W so that the pair (W, W) 
is exchangeable. This pair is then used cleverly to estimate differences in expectations of 
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test functions with respect to W and some standard distribution. Since the symmetries 
in Theorem [7| are continuous rather than discrete, it is possible to make an "infinitesimal" 
change in W by making a small change, scaling appropriately, and then taking a limit. 

Proof of Theorem^ Suppose that X is spherically symmetric and isotropic. By the spherical 
symmetry, we may assume 9 — e\, so W — X\. 

To prove the theorem, it suffices to bound |E/(W) — Ef(Z)\, where Z is a standard 
normal random variable and / : R — > R is smooth with compact support. (The proof does 
not require this much regularity of /, but it produces no loss in generality.) Recall that the 
standard normal distribution is characterized by the identity 

E[g'{Z) - Zg(Z)} = 

for all sufficiently regular g, and that given a test function /, 

g(t) = UJ(t) = f [f{s) - Ef(Z)]e-^ 2 ds 

J -co 

satisfies the differential equation 

g'(t)-tg(t) = f(t)-Ef(Z). 

In particular 

(24) E[g'(W) - Wg{W)} = Ef(W) - Ef(Z). 

To carry out the infinitesimal rotations idea described above, define a family of random 
variables {VK £ }, for e G (0, |) as follows. Let A e be the n x n orthogonal matrix 



ln-2- 



-e VT^T 2 

Now let U be a random n x n orthogonal matrix, chosen independently of X according to 
Haar measure; U T A £ U is a rotation in a random two-dimensional subspace through an angle 
sin~ 1 (5). Define 

W £ = {{U T A £ U)X,e 1 ). 

By the rotational invariance of X, (W, W e ) is an exchangeable pair for each e. 
The following facts about the joint distribution of (W, W £ ) will be needed: 

(25) E[W- W £ \W] = (l + 0(e 2 ))— W, 



n 
2e 2 



(26) E[(W - W £ ) 2 \W] = (1 + 0{e))^E[Xl\W} + 0{e A )W 2 , 

(27) E\W- W £ \ 3 = 0(e 3 ). 

Here and throughout this proof the O notation refers to asymptotic behavior as e — > 0, with 
deterministic implied constants (that may depend on n, f, or the distribution of X). The 
proof of (J23j) is given below. The proofs of (j2T?j) and (|2~7jl are similar; analogous estimates 
are proved in detail in [2*T] . 

First observe that by exchangeability, 

E((W-W e )[g(W)+g(W £ )))=0, 
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because the expression is antisymmetric in W and W £ . (In the language of Section this is 
essentially the observation that ETa = 0, where a has been chosen as in the remark at the 
end of Section El) Now by Taylor's theorem, 

(W - W £ )[g(W) + g(W £ )} = 2(W - W £ )g(W) + (W - W £ )[g(W £ ) - g(W)} 

= 2(W - W £ )g(W) -(W- W £ ) 2 g'(W) + R, 

where 

\R\ < l\W - WeWWoo <\W- W £ \ 3 \\f" 



OO 



by ©. 

By (EHD, ©, (EZD, and (HQ), 

71 

= —E[2g(W)(W - W £ ) - 2(W - W £ ) 2 g'(W)] + 0(e) 

= — E [E[2g(W)(W - W £ ) - 2(W - W £ ) 2 g'(W)\W]] + 0(e) 
= E [Wg(W) - E[X 2 2 \W]g'(W)} + 0(e), 
and so, letting e — > 0, 

E[Wg(W)]=E[E[X 2 \W]g'(W)]. 

Therefore 

(28) E[f(W) - f(Z)] = E[g'(W) - Wg(W)} = E [(l - E[X 2 2 | W])</(W)] 

for any smooth, bounded /. In particular, 

(W(£(W0, 7M ) = sup\Ef(W)-Ef(Z)\ < 4E|1 -E[X 2 2 |iy]|, 
/ 

where the supremum may be taken over smooth, compactly supported / with ||/||oo < 1, so 
that ||</||oo < 4 by (fT2"|) . This proves the first estimate. 

To prove the second and third estimates, observe that by spherical symmetry, 

E[X 2 \X,] = -1_ (EfllXH 2 !^] -X 2 ) , 



and therefore 



E|l-E[X 2 2 |Xi]| < -^(Eln-EfUXH 2 !^] | + E\X 2 - l|) 



n 




1 




1 




n 




1 




1 




n 




1 



< 7 JVar(||X|| 2 ) + 2 



by the Cauchy-Schwarz inequality and the isotropicity of X 
Finally, to prove (J2HJ), first observe that 

A £ = I n + £ J- (l + 0(e 2 ))^-I 2 



Jn-2, 
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where n _2 denotes the (n — 2) x (n — 2) zero matrix and 

Denote by K the 2 x n matrix consisting of the first two rows of the random orthogonal 
matrix U. Then 

-2 

{U T A £ U)X = X + K T ' 



eJ-{l + 0(e 2 ))-I 2 



KX, 



and so 

(29) W - W £ = -e ((K T JK)X, e x ) + (l + 0(e 2 )) ^ ((K T K)X, e x ) . 

Now if Uij denote the entries of U, then by expanding in components, 

n 

E [((K T JK)X, ei )\X] = ^^(MuiiM - u 21 u lt ) 

i=2 

and 

n 

E [((K T K)X, ei) \X]=Y1 X * E (uuu u + u 21 u 2i ). 

i=l 

Computing these expectations is not difficult because the distribution of U is unchanged by 
multiplying any row or column by —1, and any row or column of U is distributed uniformly 
on S n ~ 1 . Therefore ^ 

EuijUki = 5 ik 5ji—, 
n 

and so 

E [((K T JK)X,e 1 )\X] = 

and 

E [((K T K)X, ei )\ X] = -W. 

Putting these together proves (|25p. 

The proofs of ()26|) and (|27jl follow similarly from ()29p. One needs in addition that E|Xj| 3 < 
oo and that E[Xi|Xi] = and E[X 4 2 |Xi] = E[X||Xi] for % ^ 1. One also needs values of 
fourth order moments of the entries of U, which can be found, e.g., in |2j. 

□ 



The proof of Theorem is essentially the same as the proof of Theorem [7| the only 
difference is that A e is conjugated by a random unitary matrix instead of a random orthogonal 
matrix. The relevant mixed moments of entries of a random unitary matrix can also be found 
in0. 

Proof of Corollary^ If X is uniformly distributed on the sphere of radius y/n, then 

Var (||X||1) = 0, 

and the result follows immediately from Theorem If X is uniformly distributed on the 
ball of radius \Jn + 2, then it is easy to show by integration in polar coordinates that 
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Var < 4, from which the result follows. In both cases the constants, although 

not the order in n, can be improved by working directly from the first estimate in Theorem 
□ □ 



Proof of Corollary \ Hi By applying the Poincare inequality © to the function f(x) = 
one obtains 

Var(||X©<^, 

and so by Theorem 



By testing || on a linear functional /, one obtains that Ai < 1 when X is isotropic, and 
therefore 

d TV (L(W e ),-f R ) < 



^(v^-l) 

Since dxv < 2 always, this gives a result only for n > 25, for which the stated estimate now 
follows. □ 
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